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ABSTRACT 


A  planar  time -optimal  formulation  is  used  to  study 
the  terminal  phase  of  Interception  above  the  atmosphere 
in  a  uniform  gravitational  field.  The  dependence  on 
initial  conditions  of  the  optimal  steering  angle  is 
examined  for  initial  relative  velocities  ranging  from 
5,000  ft/sec  to  50,000  ft/sec  and  initial  distances 
up  to  500,000  ft.  Results  are  presented  in  graphical 
form  for  two  typical  rockets  showing:  (l)  the  range 
dependence  of  terminal  error  sensitivities  to  errors  in 
measurements  of  initial  conditions  taken  along  the  initial 
rectilinear  coasting  path  and  (2)  the  variation  in  these 
terminal  error  coefficients  along  optimal  interception 
paths.  Some  interceptions  with  proportional  control 
systems  were  made  and  compared  with  the  optimal  paths. 


This  report  has  been  reviewed  and  is  approved. 
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SECTION  1 


INTRODUCTION 

We  examine  the  terminal  phase  of  an  interception 
problem  in  terms  of  a  minimum  time  formulation.  The 
problem  is  based  on  a  simplified  model;  the  target  is 
assumed  to  be  non -maneuvering  and  moving  In  a  plane 
tender  the  influence  of  a  'uniform  gravitational  field. 

The  entire  interception  takes  place  above  the  sensible 
atmosphere.  The  minimum  time  interception  steering 
program  is  a  constant  angle  in  inertial  space.  The 
simple  form  of  the  optimal  control  is  appealing;  there 
may  be  situations  where  time  optimality  is  unimportant, 
yet  the  form  of  the  control  law  may  prove  to  be  a  signi¬ 
ficant  factor  in  the  design  of  the  flight  control  and 
guidance  systems.  The  terminal  phase  of  an  interception 
may  take  place  in  a  very  short  span  of  time.  For 
example ,  only  about  two  seconds  elapse  between  the  time 
a  target  is  acquired  and  the  time  a  target  is  passed  if 
the  Initial  range  is  20  miles  and  the  relative  speed  Is 
50,000  ft/sec.  In  this  study  we  do  not  try  to  map  out 
the  regions  of  accessibility  for  different  interceptor 
rocket  parameters.  We  do  present  some  typical  trajec¬ 
tories  so  that  the  nature  of  the  optimal  paths  can  be 
seen  in  some  specific  instances.  Our  goal  in  the  work 
was  to  see  how  the  optimal  control  could  be  used  In  an 
operational  sense,  to  point  out  the  critical  areas  from 
the  standpoint  of  feasibility,  and  to  provide  an 
analysis  and  a  computer  program  for  simulated  inter¬ 
ceptions.  Relative  velocities  ranging  from  5000  ft/sec. 
to  50 * 000  feet,  and  initial  (no  control)  miss  distances 
up  to  43,000  feet  are  used  In  the  illustrative  examples. 


Manuscript  released  by  the  authors  July  1964  for  publi¬ 
cation  as  an  RTD  Technical  Documentary  Report. 
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The  error  coefficients,  i.e.,  the  partial  derivatives 
of  the  final  miss  distance  with  respect  to  initial  conditions, 
form  an  essential  part  of  the  study.  These  coefficients 
depend  on  all  the  Initial  conditions.  The  computation 
of  the  optimal  control  depends  on  the  use  of  measured 
values  of  the  relative  range  and  velocity  between  the 
target  and  interceptor.  These  measurements  may  be 
accomplished  by  ground  based,  or  missile  based  equipment. 

The  control  computer  uses  these  measurements  of  relative 
position  and  velocity  to  calculate  the  optimal  steering 
program.  Errors  in  measuring  position  and  velocity 
generate  errors  in  the  calculated  steering  program 
which  in  turn  cause  errors  in  the  path  and  final  miss 
distance.  The  interceptor  must  have  a  means  for  guiding 
the  direction  of  the  rocket  thrust.  Errors  in  carrying 
out  the  commanded  steering  program  also  cause  errors  in 
the  end  conditions.  The  sensitivity  of  the  entire 
system  operation  to  random  errors  occuring  at  various 
points  in  the  control  loop  must  be  determined  so  that 
steps  may  be  taken  to  minimize  the  expected  miss  distance, 
variance  of  the  terminal  error,  or  some  other  performance 
criterion.  We  will  not  consider  the  problem  of  optimizing 
statistical  properties  of  the  end  conditions.  We  will 
outline  a  method  for  estimating  the  statistics  of  the 
terminal  errors. 

The  control  for  a  time  optimal  interception,  based 
on  a  deterministic  formulation,  ignites  the  rocket  as 
soon  as  the  target  is  sighted.  A  control  policy  based 
on  stochastic  considerations  might  not  do  this.  The 
initial  estimates  of  relative  position  and  velocity  may 
be  poor  and  might  be  expected  to  improve  as  more  obser¬ 
vations  are  made  along  the  coasting  path.  That  is,  it 
may  be  better  to  continue  coasting  towards  the  target 
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making  measurements  and  waiting  for  an  opportune  time  to 
ignite  the  rocket.  It  still  remains  to  specify  what  is 
meant  by  an  opportune  time.  Suppose  that  once  the 
optimal  steering  angle  is  determined  the  rocket  direc¬ 
tion  is  held  fixed  in  space  by  a  stabilization  system 
and  further  corrections  are  not  allowed;  this  is  an 
open-loop  terminal  guidance  scheme.  The  best  time  to 
fire  is  the  time  when  the  sensitivity  to  errors  in  the 
initial  conditions  is  least  provided  that  the  intercept 
can  still  be  made.  The  definition  of  least  sensitivity 
to  errors  is  open  to  interpretation,  but  this  will  not 
be  pursued  in  any  greater  detail  at  this  point.  Correc¬ 
tions  to  the  steering  angle  made  during  powered  flight 
require  a  closed-loop  terminal  guidance  system.  The 
sensitivity  to  °rrors  depends  not  only  on  the  sensitivity 
coefficients  taken  at  the  initial  time  but  on  the  inte¬ 
grated  effect  ox’  errors  made  during  powered  flight. 
Several  sample  sets  of  open-loop  error  coefficients  are 
computed  along  coasting  paths,  as  well  as  closed-loop 
coefficients  computed  along  optimal  trajectories.  These 
error  coefficients  yield  the  expected  terminal  errors 
for  both  open-loop  systems  and  closed-loop  flight  control 
systems  for  the  sample  optimal  paths.  The  closed-loop 
analysis  requires  knowledge  of  the  noise  characteristics 
of  the  sensors  and  other  components  of  the  guidance 
systems . 
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SECTION  2 


FORMULATION  OF  THE  PROBLEM 
I.  Time  Optimal  Interception. 

me  ioXfnpj.i.i  xeo.  t^quutionS  oi  hiutiuri  axfe  c  i-  i o // r i  ut iOW 

In  canonical  form 


v 


x 


m,/p  -  t 


m 


»c/P  - 


cos 

sin 


a 

a 


y 


Figure  1.  Coordinate  System 

We  want  to  bring  the  state  variables  x  and  y  to  zero 
simultaneously  in  the  least  time.  There  are  no  specifi¬ 
cations  placed  on  the  terminal  values  of  v  and  v  . 

x  y 

In  the  equations  above,  cQ  is  the  rocket  character¬ 
istic  exhaust  velocity,  m0  is  the  initial  system  mass, 
p  is  the  propellant  mass  flow  rate,  end  a  Is  phe  direc¬ 
tion  of  the  rocket  thrust  measured  from  the  positive  x 
axis .  The  coordinate  system  is  shown  in  Fig,  1,  The  target 
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Is  assumed  to  be  falling  freely  under  a  uniform  gravitational 
field.  The  origin  of  the  coordinate  system  is  placed  at 
the  target.  The  entire  Interception  Is  assumed  to  take 
place  above  the  atmosphere. 


The  optimal  control  is  easily  obtained  from  the 


maximum  principle  (l,  2).  Introduce  adjoint  variables 
p^(i  *  u.,.4)  and  form  the  Hamiltonian 
4 


H  -  S  P1i1  -  p 


+  P2y  +  p3vx  +  Pi)Yy  +  po 


H  "  Plvx  +  P;?vy  +  p3n7p^T  cosa  +  P.W 


—  cos  a  + 


t  sln  a  +  PG 


Maximizing  H  over  all  a  yields  the  result 


Tan  a  - 


P4O 

p3i-i 


The  adjoint  variables  p±  satisfy  the  equations 


Px(t)  =  P2(t)  s  1 
P3(t)  =  -  p1(o) 
p4(t)  =  -  p2(0) 


Therefore 

Tan  ct° 


-P-,  (0) 

-  = 


constant 


The  optimal  angle,  a0  ,  and  the  optimal  time  t° 
for  each  initial  condition  can  be  found  by  solving  a  pair 
of  transcendental  equations.  These  equations  are  obtained 
from  the  integration  of  the  equations  of  motion: 


x(t0>O  4J]coS  a°  +  vx(0)t°  +  *(0> 

(^1-^2).  ^ejain  o°  +  vy(0)t°  +  y(0) 
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These  equations  may  be  solved  by  squaring,  eliminating 


a0 


and  finding  t 


The  relations  between  t 


o 


a 


and  the  initial  conditions  are  then  used  to  find 


a 


o 


A  graphical  solution  developed  by  Faulkner  is 
described  in  Reference  2.  The  method  used  in  this  study 
is  described  in  the  Appendix. 

2.  Closed-Loop  Systems  with  Random  Noise. 

The  solution  of  the  equations  developed  in  Part  1 
of  this  section  gives  the  optimal  steering  angle  as  a 
function  of  the  initial  conditions.  The  closed-loop 
flight  control  system  based  on  this  steering  law  will 
have  to  operate  in  the  presence  of  noise.  We  distinguish 
between  the  actual  state  of  the  system,  the  measured 
state,  and  the  best  estimate  of  the  state;  the  control 
is  computed  on  the  basis  of  the  best  estimate  of  the 
state.  If  no  filter  is  used,  the  control  is  computed 
directly  from  the  measured  state  variables.  The  actual 
trajectory  evolves  from  the  initial  state,  x(0)  ,  in 
response  to  the  applied  control.  With  no  feedback  the 
trajectory  x(t)  would  miss  the  target  by  an  amount 
calculable  from  the  equations  of  motion  and  the  error 
in  control  application.  The  information  on  measurement 
errors,  disturbances,  and  control  actuation  errors  is 
assumed  to  be  statistical,  i.e.,  we  may  not  know  in 
advance  what  the  errors  will  be  but  we  may  know  something 
about  the  probability  density  distribution  of  the  errors. 
Let  the  difference  between  the  nominal  trajectory  x°(t) 
and  the  actual  trajectory  x(t)  be  called  fix  i.e.: 

.o, 


x(t)  =  X  (t)  +  fix 
motion  are: 


The  equations  governing  the 


*i(t)  -  fA(x  ,  a) 


(1) 
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where  x  is  an  n-vector  and  a  is  a  scalar  control. 


We  assume,  at  least  for  the  present,  tha 
ent table  with  respect  to  x  and  a.  We 
that  the  optimal  control  a°(x)  can  be 
everywhere  with  respect  to  x.  We  shall 
optimal  nominal  paths  to  those  which  car, 
a  field  and  eliminate  a  large  class  of  p 
bang-bang  minimum  time  controllers .  It 


t  f  is  differ- 
will  also  assume 
different iated 
restrict  the 
be  embedded  in 
roblerns,  e.g., 
is  possible  to 


analyze  the  bang-bang  problems  separately  and  they  will 


not  be  considered  here  at  all.  Kalman  has  studied  prob¬ 
lems  of  optimal  filtering  arid  control  for  linear  systems 
with  quadratic  performance  criteria  (3,^).  Our  analysis 
bears  a  resemblence  to  some  of  that  work,  but  the  motiva¬ 
tion  here  is  different.  We  want  to  estimate  the  statis¬ 
tical  properties  of  the  terminal  error,  but  we  are  not 
going  to  synthesize  an  optimal  filter. 

We  take  up  the  case  of  no  filtering.  We  want  to 

look  at  the  difference  between  the  actual  path  and  the 

nominal  path  and  eliminate  the  observed  variables.  The 

observed  values  are  x  =  x(actual)  +  N  -  x°  +  6x  +  N 

where  N  Is  additive  noise.  The  control  applied  is 

a  =  a  +  6a  where  6a  -  -r —  ~(6x  +  N)  +  a.  The  error 

dx 

in  applying  the  control  is  a.  The  sketch  below  Illus¬ 
trates  the  general  idea. 


Figure  2,  Relation  between  Actual  Path, 
Nominal  Path  and  Observed  Path 
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All  partial  derivatives  are  evaluated  along  the  nominal  path. 
The  jxo  will  not  be  used  unless  there  is  a  possibility  of 
ambiguity.  it  follows  from  the  equations  of  motion  that 


ox 


Af 

17 


If 

^  Cl 


6  a 


(2) 


Substituting  for  6a 

we  finH  that 

If 

If 

~~(6x  +  N)  +  7] 

L  dx 

(3) 

DX  *  §; 

OX  +  3— 

da 

Rearranging  terms  we  get: 


6  x  = 


df  6a° 

TS  17j 


5x 


If  1a° 
+  la  lx 


N 


(4) 


The  response  of  the  system  to  initial  errors,  measurement 
errors,  and  actuation  errors  is  evaluated  by  solving 
Equal ion  (4). 


The  changes  in  control,  6a  ,  depend  linearly  on 


the  actual,  disturbances  in  the  state,  i.e.,  6a  = 


The  term 


Jq° 

Sx 


oau 

3x " 


6x 


is  just  a  time-varying  gain.  If  we 


x 


o 


dac 


have  already  synthesized  an  optimal  controller  and 
1s  available  we  can  examine  the  behavior  of  closed-loop 
optimal  systems.  The  analysis,  however,  does  not  depend 
on  having  an  optimal  control.  It  requires  only  a  nominal 
path  and  a  linearized  feedback  equation. 


This  means  that  Equation  (4)  can  be  used  to  study 
any  linearized  feedback  system  with  time-varying  gains. 
Furthermore,  there  is  the  possibility  of  selecting  these 
time-varying  gains  to  optimize  a  performance  criterion 
associated  with  the  closed-loop  errors.  This  is  the 
approach  Kalman  takes  in  the  previously  mentioned  work 
on  optimal  linear  feedback  systems  with  quadratic  criteria. 
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Two  examples  are  given  Vo  sh-  w  now  ‘ills  method  e a n 
be  applied. 

Example  1.  Linear  Servomechanism  (  regulator*}  . 

The  system  Is  described  by  the  following  equation 

x  =  Ax  +  Bu 

w  r* 

x  is  an  n~ vector* 

A  is  n  x  n  matrix 
B  is  n  x  r  matrix 
u  is  an  r -vector 

The  commanded  control  is  a  linear  function  of  the  observed 
state  variables  u  =  kx  =  k(x°  +  6x  +  N).  The  control 
actually  applied,  u,  is  assumed  to  be  the  commanded 
control,  u,  modified  by  additive  noise,  i.e.  u  =  u  +  u, 
where  u  is  the  error  in  the  control.  The  vector  N  is 
the  noise  added  to  the  state  variables.  The  actual  errors 
In  the  state  variables  are  obtained  by  substitution. 

6x  =  +  B 

a)  If  u  s  0,  i.e.  perfect  control,  and  N  =  0,  no 
noise,  then  6x(t)  Is  determined  by  6x(0)  only, 

b)  If  6x(0)  -  0,  then  6x(t)  Is  determined  by  the 
noises  N  and  u.  The  linear  case  with 
Gaussian  distributions  for  the  random  errors 
can  be  solved  completely  by  well-known  methods. 

Example  2.  A  Time -Optimal  Steering  Problem. 

Suppose  we  desire  to  steer  a  vehicle  with  constant 
speed,  v,  from  one  location  to  another  in  least  time. 

The  equations  of  motion  are 

x  =  v  cos  a  (6a) 

y  =  -v  sin  a  (6b) 

where  a  is  the  angle  between  the  velocity  vector  and 
the  x-axls  measured  positively  counterclockwise. 
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+  B 


c)x 


N  +  Bu 


(5) 


3 


Suppose  th‘-  ‘.arret  j  s  a  I  r  ■  i  of  rad  1  u r>  '■<  -entered 

a.t  the  origin  as  shown  in  Figure  3  below  . 


Figure  3*  Time  Optimal  Steering  Example 

It  is  easy  to  show  that  tne  time  optimal  steering 
angle  for  this  problem  directs  the  velocity  vector  towai’ds 
the  target  at  each  instant  of  time.  The  isochrones  are  circles 
centered  at  the  origin.  The  application  of  Equation  (4)  yields 


The  equation  of  interest  is  (7b).  We  are  free  to 
choose  the  orientation  of  the  coordinate  system  axes. 
The  nominal  initial  conditions  can  be  taken  as: 
x(0)  =  xq,  y(0)  =  0  and  the  nominal  optimal  path  will 
be  y°(t)  £  0  . 


ba° 


Along  the  nominal  path 
.-1  y 


lyC 


sin 


ba° 

6x‘ 


y 

y° 


*1 


=  0 


and 


=  -  - 4  for  xd»yd  where  R2=x2+y2 


R 


The  equation  for  5^  becomes 
^  (6y+N2)  +  VO 
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x  (O)-vt 


vw2  +  va 

X  (  fcf}"-  vt 


The  equation  for  the  closed-loop  optimal  control  has  a 
singularity  at  the  nominal  interception  time;  t  =  • 

The  net  effect  is  to  create  a  feedback  gain  inversely 
proportional  to  R. 

If  the  initial  error  6y(o)  is  known  and  if  the 
noise  Ng  and  actuation  error  are  known  to  be  statistically 
independent  then  the  statistical  behavior  of  6y (t )  can 
he  found  from  the  Fokker-Planck  equation  (5)*  In  the 
case  of  Gaussian  noise,  the  solutions  are  easily  obtained 
by  more  direct  methods.  For  example,  the  solution  of  the 
differential  equation  can  be  written  explicitly  as 


fiy(t)=&y(o 


V")dT 


/xl 

(O)-vt^ 

i  i 

r  xt°) 

V  x(o;  j 

J 

x(0) 

/■>  (.x(O)-vt)' 

o 

a(x  )dT 


Suppose  =  K(t)  where  K(t)  is  prescribed  in 

advance  Instead  of  the  optimal  feedback.  The  equations 
then  are 

6$  =  vK(t)  (6y+Ng)  +  va 
If  K(t)  is  constant,  then 

6f  =  VK(6y+N2)  +  va 


and  the  steering  problem  now  is  a  special  case  of  the 
servomechanism  problem  used  in  Example  1  with  n  *  1. 
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This  first  order  differential  equation  represent;.: 
a  linear  transformation  of  Gaussian  processes.  The 
random  function  5y(t)  is  the  sum  of  two  lirearly  trans¬ 
formed  Gaussian  processes  and  has  a  Gaussian  distribution. 

Suppose,  for  example,  that  N2  =  C.  Then,  return¬ 
ing  to  Example  2, 


o 

t 


where  kQ^Ti,T2^  “  E^(a(T-,)  -  a)(a(r2)  -  a)|  is  the 
autocorrelation  function  for  a  (6),  The  same  result 
can  be  obtained  via  the  weighting  function -power 
spectral  density  approach  (7). 

The  analysis  discussed  in  the  preceding  paragraph 
can  be  applied  in  a  relatively  simple  way  to  improve 
non-optimal  feedback  about  an  optimal  path.  Suppose 
that  we  use  a  constant  gain  feedback  (K  <  0)  based  on 
the  displacement  from  the  nominal  path.  If  the  actuator 
errors  are  negligible  compared  to  the  measurement  errors, 
we  may  write:  by  a  vK6y  +■  vKN. 

Suppose  also  that  the  initial  errors  and  measure¬ 
ment  noise  have  Gaussian  distributions  with  zero  means. 
Let  the  variance  of  the  measurement  noise  be  q,  and 
the  noise  be  "white"  such  that  k^(xp T2^  ~  ^ ^ Ti  *  fg) 
and  5  Is  the  Dirac  Delta  function;  the  variance  of  the 
initial  errors  be  a  (0).  The  variance  at  the  terminal 
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time  t°  is  given  by 


a 


2 

6y( t°) 


.=  a 


2  2vKt° 
6y(0)e 


3^(i 


e2vKt°) 


and  6y(t°)  has  a  Gaussian  distribution  with  zero  mean 

p 

and  variance,  g  ,  c,  .  It  is  possible  to  find  a  value 

^  y  °  *  9  ^  o 

of  gain  K  that  minimizes  c ,  t  o\  :  no*. e  *  r.a^ 

J  S'  ■*  "  i  f  !  A  ry  /■',  v  /■  »  t  «. 

^y  \  Or.  J 

<  0  at  K  -  0  ,  and  of  ,-^o\  — ^  +0D  as  K  — u  -hoG  .  We 

<jy  !o  j 


will  not  carry  these  simple  examples  any  further.  They 
serve  the  purpose  cf  showing  that  the  error  coefficients 
are  desirable  objects  to  obtain.  They  play  an  important 
role  in  evaluating  the  behavior  of  a  control  system  in  the 
presence  of  random  noise.  In  the  next  section,  we  will 
examine  the  dependence  of  the  error  coefficients  on 
initial  conditions. 


SECTION  3 


COMPUTATIONAL  STUDIES 

The  discussion  in  Section  2  was  directed  at  the 
analytical  formulation  of  the  interception  problem.  In 
this  section,  we  examine  the  sensitivities  of  the  inter¬ 
ception  paths  to  errors.  Figure  4  shows  a  typical  inter¬ 
ception  configuration.  The  parameters  of  interest  are 
the  initial  miS3  distance,  pQ  =  R(0)  sin  0,  the  initial 
rarge  R(0),  and  the  initial  velocity  V(0).  The  thrust 
level  is  not  adjustable;  the  optimal  steering  angle  can 
be  found  by  the  method  described  in  the  appendix. 


Table  I  contains  a  partial  listing  of  the  cases  studied. 

The  miss  distances  or  distances  of  closest  approach,  p, 

are  given  for  five  values  of  initial  range  and  four  values 

of  initial  velocity.  The  results  are  grouped  according 

to  the  initial  range;  each  group  is  ordered  according  to 

the  initial  miss  distance  p  . 

o 

Two  rockets  were  used  In  this  study  and  their 
characteristics  are  listed  below.  These  values  are 
typical  of  the  kina  of  rockets  being  considered  for 
use  In  missile  and  satellite  interception  studies. 
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-rz?  -avi 


Rocket  #1 


cQ  =  10,000  ft/sec-  effective  exhaust  velocity 
=  .15  mQ,  3/mQ=  .1/sec,  tb  =  8.5  sec 
a  (0)  =  30  g  (approx. ) 
a  (tfc)  =  200  g  (approx.) 

AVideal  =  18^00  ft/sec 


Rocket  #2 

cQ  =  10,000  ft/sec 

rrij.  =  .15  mQ,  0/ m  =  .05/sec,  tfe  =  17.0  sec 
a  (0)  =  15  g  (approx.) 
a  (tb)  =  100  g  (approx.) 

AVideal  =  18^900  ft/sec 

These  rockets  have  the  same  mass  ratio  (m£./m0)  and  the 
same  ^ .  They  differ  In  the  burning  rate  of  the 

propellant  and  in  the  acceleration  time  history. 

Two  sets  of  trajectories  are  shown  in  Figures  6 
and  7  .  (Note  that  the  x-y  scales  are  not  the  same.) 
Each  set  of  figures  shows  the  effect  of  initial  speed 
on  the  interception  paths;  Figure  6  is  for  Rocket  #1, 
Figure  7,  for  Rocket  #2.  It  may  help  to  refer  back  to 
Figure  1  where  the  coordinate  system  is  defined. 

The  trajectories  shown  in  Figures  6  and  7  represent 
simulated  interceptions  with  the  rocket  ignited  at  the 
point  of  acquisition  of  the  target. 

Imagine  a  situation  where  there  is  a  possibility 
of  delaying  the  Ignition  of  the  rocket  to  Improve 
the  relative  position  and  velocity  measurements. 


Figure  8  shows  the  family  of  interception  paths  for 
Rocket  #1  emanating  from  the  coasting  path  starting  at 
R(0).  The  target  is  sighted  at  a  range  of  500,000  feet 
and  the  interceptor  approaches  the  target  on  the  coast¬ 
ing  path  at  a  speed  of  10,000  ft/sec.  The  miss  distance 
p  in  the  absence  of  terminal  control  is  4,363  feet. 

Each  path  in  this  family  is  obtained  by  solving 
for  the  optimal  steering  angle  and  computing  the  optimal 
trajectory.  The  optimal  steering  angles,  a°,  are 
shown  in  Figure  8  next  to  the  starting  points.  The 
optimal  steering  program  to  the  target  from  any  start¬ 
ing  point  is  constant  and  is  measured  from  the  positive 
x-axis  as  shown  in  Figure  1.  This  angle  depends  on  the 
initial  conditions  and  it  changes  as  the  starting  points 
vary  along  the  initial  coasting  path.  The  optimal  steer¬ 
ing  angle  changes  by  about  only  7°  as  the  initial  range 
decreases  from  500,000  feet  to  about  100,000  feet.  The 
rate  of  change  of  steering  angle  Increases  sharply 
thereafter.  The  effect  of  a  time  delay  in  igniting  the 
rocket  can  be  determined  by  an  examination  of  the  rate 
of  change  of  a°  along  the  coasting  path. 

Open-loop  miss  sensitivities  have  been  found 
for  each  of  the  starting  points  shown  in  the  figure. 

These  partial  derivatives  represent  the  sensitivity 
of  terminal  errors  to  errors  in  initial  measurements 
and  to  t-rrors  in  steering.  These  are  shown  next  in 
Figure  9.  The  sensitivity  of  the  error  in  y(T)  at 
the  distance  of  closest  approach  to  errors  in  a° 
is  shown  in  Figure  9*  The  sensitivity  of  y(T)  to 
errors  in  a°  decreases  continually  as  the  target 
is  approached  until  the  range  is  very  small.  A  sharp 
Increase  in  steering  angle  sensitivity  with  decreas¬ 
ing  range  is  shown  for  the  other  initial  conditions  in 


16 


the  remaining  figures  of  tnis  set.  This  is  not  unex¬ 
pected;  the  a°  listings  or.  Pig.  8  show  a  rapid  change  in 
the  initial  steering  angle  as  the  range  approaches  40,000 
ft.  The  paths  change  from  head-on  interceptions  to  beam 
approach  interceptions  in  this  region  as  the  interceptor 
coasts  by  the  target.  The  initial  relative  velocity  of 
10,000  ft/sec  is  low  enough  to  permit  the  interceptor  to 
reverse  its  direction  of  motion.  Two  paths  shown  in  Fig. 8 
pass  the  target  and  reverse  direction.  The  discontinuous 
changes  in  the  error  sensitivities  (shown  in  Fig.  9)  occur 
as  the  interceptor  passes  a  beam  of  the  target. 

It  is  possible  to  compute  some  of  the  partial  deriva¬ 
tives  directly  from  the  integrals  of  the  equations  of 
motion,  e.g.. 


c>y(o) 


These  derivatives  are  evaluated  at  the  nominal  optimal 
o  £qo  aao 

time,  t  .  The  derivatives  of  the  form  j  5£(6t  *  etc • > 

are  more  difficult  to  obtain  in  closed  form  as  the  optimal 

angle  is  given  implicitly. 

Trajectories  with  Initial  relative  speeds  of  25,000 

ft/sec  and  50,000  1't/sec  are  shown  in  Figures  10  and  12 

with  Rocket  #1  and  in  Figures  14  and  16  with  Rocket  #2. 

Error  coefficients  for  these  cases  are  given  In  Figures 

11,  13,  15,  and  17. 
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The  principal  error  coefficients 


dy(T)  ^  ayj  T)  ?  dy(T) 
3x(0)  ’  5xf0)  ’  dy(O) 


are  presentee  .in  hugurcs  13  and  1 ')  ,  or  toe  trs.ifec  uor.i  cs 
of  Figures  10,  12,  l4,  and  10.  These  are  evaluated  at  T  , 
the  time  cf  closest  approach  to  the  target,  with  perturbed 
thrust  direct  Ion  angle  res  iting  from  perturb  at inns  In  the 
initial  conditions.  The  error  coefficients 


dy(T)  f  by(7) 

Sx(0)  bx(o) 


increase  in  absolute  magnitude  as  the  initial  range 

bv^) 

decreases;  the  coefficient  decreases  in  aosolute 

dy(  0) 

magnitude  as  the  range  decreases.  Trie  location  of  the 
best  place  to  ignite  the  rocket  from  the  standpoint  of 
open-loop  terminal  error  dispersion  can  be  found  by  using 
the  noise  characteristics  of  the  sensors. 

The  analysis  of  feedback  systems  requires  the 
partial  derivatives 


ba°  bo°  da0 

bx  ’  bx  by 


and 


da 


along  the  optimal  path.  The  latter  three  derivatives  are 

presented  in  Figures  20  and  21  for  theosainple  cases 

discussed  earlier.  The  derivative  -  is  very  small 

ox 

except  at  the  very  end  of  the  flight  and  is  rot  plotted 
here.  The  derivatives  are  plotted  as  functions  of  R  , 
range.  The  differential  equations  require  these  as 
functions  of  time  and  the  computer  program  provides  both 
forms.  The  derivatives  generated  by  the  computer  program 
can  be  used  as  indicated  in  the  following  paragraphs. 
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If  the  equations  for  error  propagation  in  an 
optimally -controlled  system  are  expanded,  we  obtain 

<$x  =  6v„  i,  i. 


I 

6vx=  -a(t)  sin  a°[^||p5x1 

i=l  1 

«5y  =  *vy  4 

6vy=  a(t)  cos  a0  J" 


Il^i]  '  [a(t)  sln  a°]  s 


da 

^Xi 


i=l 

4 


i=l 


■  +  [a(t)  cos  a°]  “ 
i=l 


x3  =  6y,  x4 


is  known  from  the  rocket  speci- 
o 

> 


a 


is  a  constant  along 


where  a(t)  =  t,  =  4*.  xsj.  6v 

6vy,  Ni  =  measurement  errors,  and  a  is  the  error  in 
applying  the  commanded  control. 

The  quantity  a(t) 
fications.  The  optimal  angle, 
the  nominal  path. 

The  specification  of  the  statistical  character  of 
the  noise  N.^  and  actuation  errors  a(t)  provides  the 
Information  needed  to  evaluate  the  statistical  properties 
of  the  terminal  errors.  The  function  a(t)  is  a  rapidly 

^qO  ^qO 

increasing  function  of  time.  The  derivatives  — — — , 

an0  ox  °y 

are  also  increasing  functions  of  time.  The  optimal 

steering  angle  becomes  very  sensitive  to  measurement 

errors.  This  may  be  thought  of  as  a  rapid  increase  in 

a  feedback  gain  near  the  end  of  powered  flight. 

A  time-varying  feedback  control  gain  that  takes 
Into  account  the  variation  in  measurement  error  statis¬ 
tics  along  the  flight  path  may  prove  to  be  better  than 
the  closed-loop  time  optimal  controller  from  the  stand¬ 
point  of  error  sensitivity. 

The  error  propagation  is  governed  by  a  set  of 
Inhomogeneous  time-varying  linear  differential  equations. 
The  forcing  functions  acting  on  these  equations  are  the 
measurement  errors  N^(t)  and  the  actuation  errors  a(t). 
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Large  measurement  errors  during  powered  flight  may  cancel 
the  benefits  derived  from  feedback  of  the  estimated 
position  and  velocity  errors.  If  there  were  no  measurement 
errors,  i.e.,  =  0  ,  then  the  terminal  errors  due 

to  a  (t)  can  be  reduced  by  the  use  of  feedback  of 
position  and  velocity  signals.  The  presence  of  measurement 
noise,  ,  makes  necessary  a  detailed  analysis  of 

closed-loop  guidance.  It  seems  fairly  clear  that  if 
the  noise  level  is  high  compared  to  the  initial  errors, 

6x^(0)  ,  it  may  be  desirable  to  run  the  system  open-loop. 

The  terms  closed-loop  and  open-loop  apply  to  the  use  of 
commanded  changes  in  the  steering  angle  due  to  errors 
In  position  and  velocity  along  the  nominal  path.  The 
optimal  steering  angle  in  either  case  is  to  be  controlled  or 
stabilized  in  direction  by  some  other  attitude  control  system. 

We  will  illustrate  the  use  of  results  of  the  computer 
program  by  a  simple  single  variable  open  loop  example. 

Suppose  that  the  measurements  ot  the  single  variable  y 
are  corrupted  by  additive  noise,  the  variance  .of  the 
measurement  noise  depends  linearly  on  range,  and  the 
commanded  contro]  has  no  errors  i.e.  a  s  0.  Suppose 
also  that  one  of'  the  graphs  in  Figure  18  represents  the 
sensitivity  of  terminal  errors  to  errors  in  initial 
determination  of'  position  and  that  this  curve  can  be 
satisfactorily  approximated  by  an  exponential  function. 

It  is  relatively  simple  to  find  the  best  range  for 
initiating  the  control  to  give  tne  smallest  variance  in 
the  terminal  errors.  The  best  location  occurs  either  at 
an  interior  point  or  at  one  of  tne  extreme  ends  of  the 
interval  beginning  at  the  point  or  target  acquisition 
and  ending  at  the  last  point  wnere  an  interception  can 
be  Initiated.  (This  interval  may  not  be  closed.)  We  assume 
that  the  terminal  error  depends  linearly  on  6y(0)  for 
small  initial  errors  so  that 
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CTy(T 


by  (T) 

-by(C) 

y(0)  =  y(0) 

An  interior  minimum  is  sought  by  differentiating  the 
expression  (where  y  is  used  in  place  of  y(0}) 

'  j  !  l  ^T1  i  ”] 

a";(x)  •  i  Z't'/'i'  I 

l  dy  ( O )  J 

with  respect  to  range  (x) .  In  this  example,  let 


O^(x)  =  b  +  CX 

y 


bj  c 


y  9 


x  >  0 


and 


Mil 

LSy(O) 


1  2 


=  aV2Xx 


X  <  0 


We  must  look  at  the  zero  of 


ji 
dxl 


(a^e~2Xx)  *  (b  +  cx) 


and  see  if  it  lies  in  the  interval  of  interest  and  test 
to  see  if  we  have  a  maximum  or  a  minimum.  Even  this 
simple  example  provides  some  interesting  results.  The 
derivative  of  the  terminal  variance  with  range  is 

-a2  e  '2X:xj^2X(b  +  ex)  -  c 

The  vanishing  of  the  derivative  corresponds  to  a 
maximum  sensitivity;  the  minimum  sensitivity  occurs  at 
one  of  the  endpoints.  The  minimum  location  depends 
strongly  on  the  parameter  b.  For  b  =  0,  the  derivative 
vanishes  at  x  =  ~.  The  interior  extrema  do  not  depend 
on  the  parameter  a  at  all. 

This  procedure  can  be  applied  to  multivariable  prob¬ 
lems  In  a  straightforward  way  using  co -variance  matrices 
Instead  of  the  scalar  variance  used  in  the  example. 
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SECTION  4 


PROPORTIONAL  CONTROLLERS 


Interceptions  with  proportional  control  systems 
were  simulated  on  a  digital  computer.  Three  modes  of 
control  were  studied: 


Mode  A  -  The  thrust  is  directed  along  the  line -of - 
sight  towards  the  target.  The  possibility 
of  using  a  constant  lead  angle  is  included. 


Mode  B 


The  thrust  is  directed  off  the  llne-of- 

7  v  ? 

sight  by  an  angle  proportional  to  — j-p-j — 
in  the  sense  required  to  reduce  the  rate 
of  rotation  of  the  line-of -sight . 


Mode  C  -  The  thrust  i_>  pplied  normal  to  the  line- 
of -sight  in  the  direction  needed  to  reduce 
the  rate  of  rotation  of  the  line-of -sight . 
If  the  interceptor  velocity  becomes  aligned 
with  line-of -sight,  there  is  an  option 
either  to  stop  thrusting  or  to  thrust 
along  the  line-of -sight . 


By  putting  limits  on  the  steering  angle  and  raising 
the  gain,  we  can  also  obtain  thrusting  perpendicular  to 
the  line-of -sight  from  Mode  B.  Firing  perpendicular  to 
the  line-of-sight  worked  quite  well. 

We  did  not  explore  all  the  possible  combinations  of 
initial  conditions  and  control  gains.  The  reachable  sets 
for  the  optimal  steering  were  not  mapped  in  detail.  It 
Is  not  difficult  to  determine  on  the  basis  of  physical 
considerations  which  of  the  possible  control  schemes  are 
most  desirable.  The  control  must  direct  the  relative 
velocity  vector  towards  the  origin  before  the  rocket  burns 
out.  If  the  time  of  flight  in  the  terminal  phase  is 
shorter  than  the  rocket’s  burning  time  and  if  the  trans¬ 
verse  velocity  component  is  large  compared  with  AV 
available  for  the  rocket,  then  it  is  desirable  to  fire 
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slightly  backwards.  This  lengthens  the  time  available 
for  reducing  the  transverse  component.  In  effect,  this 
is  a  trade  between  interception  time  and  terminal  error. 
Firing  along  the  line-of-sight  is  undesirable  in  this 
situation.  On  the  other  hand,  if  the  transverse  velocity 
is  low,  the  initial  miss  distance  is  small, and  the  time 
of  flight  is  long,  then  the  time  optimal  steering  angle 
is  very  close  to  the  line-of-sight.  This  can  be  seen 
in  Figure  23. 

Figure  22  shows  five  trajectories  starting  from 
identical  initial  conditions.  Path  number  1  is  the  time 
optimal  path;  thrusting  perpendicular  to  the  line-of- 
sight,  Mode  C,  gives  path  number  2.  Mode  B  produces 
path  3  with  a  miss  distance  of  24,000  feet.  The  result 
of  steering  along  the  line-of-sight.  Mode  A,  is  the  path 
marked  4,  with  a  miss  distance  of  91,000  feet.  The 
steering  angles  used  are  given  in  Figure  23. 

The  gains  used  in  these  runs  were  selected  a  priori 
and  were  not  optimized  with  respect  to  the  miss  distances. 
The  optimal  control  system  employs  a  computer  to  solve  a 
transcendental  equation;  the  application  of  this  computer 
to  the  problem  of  gain  and  lead  angle  selection  would 
undoubtedly  improve  the  performance  of  the  proportional 
controllers . 

The  optimal  controller  nevertheless  possesses 
advantages  which  are  not  matched  by  the  proportional 
controllers.  The  optimal -control  computer  will  minimise 
the  miss  distance  instead  of  interception  time  if  a  hit 
cannot  be  scored.  The  computation  of  the  optimal  control 
also  takes  into  account  the  terminal  coasting  phase  after 
the  propellant  is  exhausted.  Furthermore,  the  optimal 
control  produces  hits  for  initial  conditions  which  lie 
on  the  boundary  of  the  region  of  attainability  whereas 
the  proportional  controllers  need  not  do  this. 


SECTION  5 

SUMMARY  AND  CONCLUSIONS 


A  planar  time  optimal  formulation  is  used  to  study 
the  terminal  phase  of  an  Interception  problem.  The 
optimal  steering  program  for  interception  In  a  uniform 
gravitational  field  in  the  absence  of  aerodynamic  drag 
requires  the  thrust  direction  to  be  held  constant  in 
Inertial  space  during  the  powered  portion  of  flight. 

The  optimal  steering  angle  is  a  function  "of  the  relative 
initial  position  and  relative  velocity.  We  have  examined 
the  dependence  of  the  optimal  steering  angle  on  the  ini¬ 
tial  conditions  for  initial  velocities  ranging  from 
5,000  ft/sec  to  50,000  ft/sec.  The  results  are  presented 
in  graphical  form  showing:  (1)  the  range  dependence  of 
the  terminal  error  sensitivities ' to  error  in  initial 
position  and  relative  velocity  measurements  taken  along 
an  Initial  rectilinear  coasting  path,  and  (2)  the  varia¬ 
tion  in  these  terminal  error  sensitivities  along  an 
optimal  interception  path. 

Simulated  Interceptions  were  made  using  proportional 
controllers.  Three  modes  of  control  were  tried;  thrust¬ 
ing  perpendicular  to  the  line-of -sight  gave  good  results 
in  the  cases  studied.  The  optimal  controller  Is  predictive 
in  nature  and  takes  into  account  the  possibility  of  a 
terminal  coasting  phase  after  the  propellant  Is  exhausted. 
The  computer  also  predicts  whether  or  not  a  hit  Is  possi¬ 
ble  and  will  minimize  the  miss  distance  if  an  intercep¬ 
tion  is  not  possible.  Further  work  is  needed  to  improve 
the  performance  of  the  proportional  controllers,  e.g., 
selecting  thu  gain  and  lead  angle  as  functions  of  range, 
velocity  or  combinations  of  other  variables  before  . 
definitive  comparison  can  be  made  with  optimal  controllers. 


24 


One  of  the  results  of  our  work  is  a  flexible  digital 
computer  program  with  many  options,  which  can  simulate 
optimal  interceptions,  generate  the  error  sensitivities 
along  initial  coasting  paths  and  along  optimal  paths  and 
make  terminal  error  dispersions  studies.  A  listing  of 
the  program  and  a  set  of  cards  are  available  at  ARAP. 

The  study  of  error  sensitivities  raised  several 
important  questions  which  are  not  yet  settled.  The 
computer  program,  as  it  now  stands,  is  able  to  simulate 
closed-loop  interceptions.  Several  check  problems  were 
run  but  systematic  studies  were  not  made  because  of  the 
dearth  of  unclassified  realistic  data  on  radar  and 
sensor  noise  and  missile-borne  guidance  and  attitude 
control  systems.  Future  research  programs  in  this  area 
could  take  advantage  of  this  computer  routine  to  inves¬ 
tigate  the  behavior  of  optimally-controlled  hypervelocity 
interceptions,  with  and  without  relative  position  and 
velocity  feedback  during  the  burning  portion  of  the 
flight.  This  simulation  should  include  measurement 
errors  and  actuation  errors  based  on  experience  with 
actual  equipment  tests. 
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TABLE  1 


SUMMARY  -F  INTERCEPTION  RUNS  WITH  OPTIMAL  CONTROL 

Rocket#! 


Initial 

Range 

R(O) (ft) 

:1s s  Distance 
without 
control 

Pc.  (ft) 

Miss  Distance 
with  optimal 
control 

P  (ft) 

Initial 
Velocity 
(ft/sec ) 

v(o) 

Optimal 

Steering 

Angle 

a° 

50,000 

873 

.000 

10,000 

175.77 

It 

it 

.000 

15,000 

171.03 

M 

It 

.000 

25,000 

155.01 

tt 

If 

355.0 

50,000 

89.79 

II 

1,309 

.000 

10,000 

173.66 

ii 

ri 

.000 

15,000 

166.50 

It 

II 

.000 

25,000 

140.94 

II 

tt 

791.0 

50, 000 

90.29 

It 

1,745 

.000 

10,000 

171.54 

H 

II 

.000 

15,000 

161.91 

II 

tl 

.000 

25,000 

123.70 

II 

ll 

1,228.0 

50,000 

90.79 

II 

4,358 

.000 

10,000 

158.75 

II 

II 

.000 

15,000 

131.47 

II 

II 

2,209.0 

25,000 

89.86 

II 

tl 

3,843.0 

50, 000 

93.80 

II 

8,682 

.000 

10, 000 

136.70 

II 

If 

.000 

15,000 

13.08 

II 

It 

6,567.0 

25,000 

94.88 

II 

It 

8,179.0 

50, 000 

98.81 

125,000 

2, 182 

.000 

10, 000 

178.29 

II 

If 

.000 

15,000 

176.53 

II 

It 

.000 

25,000 

170.72 

II 

It 

.000 

50, 000 

139.36 

II 

3,272 

.000 

10, 000 

177.43 

II 

II 

.000 

15,000 

174.79 

tt 

tl 

.000 

25,000 

166.03 

II 

II 

.000 

50,000 

105.63 

II 

4,362 

.000 

10,000 

176.58 

II  | 

It 

.000 

15,000 

173.05 

It 

It 

.000 

25.000 

161.28 

Tt 

II 

936.0 

50,000 

88.70 

t! 

10,894 

.000 

10,000 

171.48 

If 

II 

.000 

15, 000 

162.64 

It 

If 

.000 

25,000 

129.56 

II 

II 

7,480.0 

50, 000 

91.70 

tt 

21,706 

.000 

10,000 

163 . 17 

II 

tl 

.000 

15,000 

145.35 

If 

6,005.0 

25,000 

89.98 

It 

18,355.0 

50,000 

96.72 
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TABLE  i  (Con’t. ) 


Initial 

Range 

R { 0 } (ft) 

Miss  Distance 
without 
control 

P0  ( ft ) 

Miss  Distance 
with  optimal 
control 

P  (ft) 

Initial 
Velocity 
( ft/sec ) 

v(o) 

j  Optimal 
Steering 
Angle 
a0 

200, 000 

3,490 

.000 

10,000 

178.89 

If 

If 

.000 

15,000 

177.94 

1 1 

It 

.000 

25,000 

174.56 

II 

.000 

50, 000 

156.96 

5,235 

.000 

io, 000 

178.34 

ir 

Tl 

.000 

15, 000 

176.92 

*  ■ 

.000 

25,000 

171.83 

?! 

ff 

.000 

50,000 

144.28 

ft 

6,  &30 

.000 

10,000 

177.78 

If 

.000 

15,000 

175.89 

H 

.000 

25,000 

169.10 

fl 

l? 

.000 

50,000 

129.52 

It 

17,431 

.000 

10,000 

174 . 47 

1! 

It 

.000 

15,000 

169.76 

tl 

It 

.000 

25,000 

152.47 

it  I 

ir 

8,067.0 

50, 000 

89.10 

ft 

34,730 

.000 

10,000 

169.03 

ir 

n 

.000 

15,000 

159.78 

tl 

tt 

.000 

25,000 

122.65 

1! 

11 

25,508.0 

50, 000 

94.12 

250, 000 

4,3(53 

.000 

10,000  1 

179.04 

n 

.000 

15,000 

178.32 

fr 

ff 

.000 

25,000 

175.86 

!l 

(1 

.000 

50,000 

162.31 

n 

6,544 

.coo 

10,000 

173.56 

it 

II 

.000 

15,000 

1  177.48 

tt 

11 

.000 

25, 000 

173.80 

tl 

tl 

.00  0 

50,000 

153.01 

tl 

8,725 

.000 

10, 000 

178.08 

tl 

ft 

.000 

15,000 

176.65 

if 

ff 

.000 

25,000 

171.73 

it 

II 

.000 

50,000 

143.03 

n 

21,789 

.000 

10,000 

175.21 

Tl 

II 

.000 

15,000 

171.64 

II 

tt 

.000 

25,000 

159.32 

II 

tt 

6, 366.0 

50, 000 

86.93 

fl 

43,412 

.000 

10,000 

170.49 

1! 

11 

.000 

15,000 

163.41 

It 

H 

.000 

25,000 

138.76 

ft 

it 

28,175.0 

50,000 

91.92 

500.000 

8,726 

.000 

10, 000 

179.28 

fl 

tt 

.000 

15,000 

178.85 

tl 

I? 

.000 

25, 000 

177.80 

ft 

f? 

.000 

50,000 

172.88 

If 

13,088 

.000 

10,000 

178.93 

II 

ff 

ft 

.000 

15,000 

178.28 

ft 

.000 

25,000 

176.70 

.000 

50, 000 

169.30 

'•'<4 


TABLE  1  (Con't.) 


Initial 

Range 

R(O) (ft ) 

Miss  Distance 
without 
control 

P0  (ft) 

Miss  Distance 
with  optimal 
control 

P  (ft) 

Initial 
Velocity 
( ft/oec ) 

v(o) 

Optimal 

Steering 

Angle 

a° 

500, 000 

17,450 

.000 

10,000 

178.57 

It 

tf 

.000 

15,000 

177.71 

II 

it 

.000 

25,000 

175.60 

ft 

tr 

.000 

50, 000 

165.71 

II 

43* 57 &  - 

.000 

10, 000 

176.42 

If 

ff 

.000 

15,000 

174.28 

ft 

If 

.000 

25,000 

168.99 

ft 

II 

.000 

50,000 

143.41 

It 

86,824 

.000 

10, 000 

172.87 

ft 

fl 

.000 

15,000 

168 . 60 

If 

It 

.000 

25,000 

157.96 

t! 

II 

.000 

50, ooo 

89.24 

Rocket#g 


50,000 

873 

.000 

10, 000 

171.71 

ft 

tf 

.000 

15,000 

161.77 

It 

ff 

.000 

25,000 

121.54 

If 

fl 

618.0 

50,000 

90.41 

If 

1,309 

.000 

10,000 

167.53 

ft 

1! 

.000 

15,000 

152.13 

ff 

tf 

273.0 

25,000 

89.08 

ff 

II 

1,055.0 

50,000 

90.91 

II 

1,745 

.000 

10,000 

163.31 

1? 

ff 

.000 

15,000 

141.71 

ft 

t! 

710.0 

25,000 

89.58 

I! 

tf 

1,491.0 

50, 000 

90.91 

ff 

4,358 

.000 

10,000 

135.85 

fl 

It 

.000 

15,000 

5.45 

It 

.  .  ....  vi-:;-.  W-v.  ur 

r  '3,32T.0  '  • 

25,000' 

92.58 

ff 

ff 

4,105.0 

50, 000 

94.41 

8,682 

.000 

10, 000 

13.83 

ft 

If 

.000 

15, 000 

10.89 

fl 

If 

7,672.0 

25, 000 

97.60 

If 

ft 

8,435.0 

50,000 

99.42 

125,000 

2,182 

.000 

10, 000 

176.60 

ff 

ft 

.000 

15,000 

172.85 

fl 

If 

.000 

25,000 

160.39 

n 

ft 

549.0 

50, 000 

89.47 

if 

3,272 

.000 

10,000 

174.90 

ft 

ff 

.000 

15,000 

169.26 

ft 

0  ff 

.000 

25,000 

149.91 

it 

ft 

1,640.0 

50,000 

89.97 

« 

4,362 

.000 

10,000 

173.19 

ft 

tt 

.000 

15,000 

165.65 

.  if 

"'V  'J  . 

tf 

.000 

25,000 

138.38 

ft 

W 

2,731.0 

50, 000 

90.47 
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TABLE  1  (Con't.) 


Initial 

Range 

R( 0) ( f 4  ) 

Ml 3G  Distance 

without 

control 

p„  (  ft ) 

u 

Miss  Distance 
with  optimal 
control 

p  (ft) 

Initial 
Velocity 
( rt/sec ) 

v(c) 

Optimal 

Steering 

Angle 

a0 

125, 000 

10,094 

.  000 

10 , 000 

162.97 

If 

If 

.000 

15,000 

142.38 

ft 

If 

}  ■  -"\  <-“N  (  "» 

H  ,  .0 

25 » Guo; 

o3 . 34 

If 

ft 

9,271.0 

50,000 

93.47 

ft 

21,706 

If 

.000 

10 , oeo 

145.39 

f? 

.  000 

15,000 

15.93 

If 

If 

14,932.0 

25,000 

93.35 

IT 

If 

20, 117.0 

50,000 

93.49 

200, 000 

3,590 

.  000 

10,000 

177.85 

ft 

.000 

15,000 

175.60 

If 

tf 

.000 

25, 000 

168.17 

I! 

It 

.000 

50,000 

124.64 

It 

5,235 

.000 

10,000 

176.78 

If 

ft 

.  000 

15,000 

173.40 

*1 

It 

.000 

25,000 

162.14 

!  1 

tf 

6, 930 

936.0 

50,000 

88.94 

11 

.  000 

10,000 

175.71 

If 

ft 

.000 

15,000 

171.20 

tl 

tf 

.000 

25,000 

155.95 

If 

ff 

2,680.0 

50,000 

89.44 

If 

17,431 

.000 

10,000 

169.31 

It 

tf 

.000 

15,000 

157.88 

If 

ff 

.000 

25, 000 

106.85 

If 

!T 

13,151.0 

50,000 

92.44 

fl 

34,730 

.000 

10,000 

158.86 

If 

ff 

.  000 

15,000 

134.95 

If 

If 

15,904.0 

25,000 

8? .  36 

It 

tf 

4,363 

30,533.0 

50,000 

97.46 

250,000 

.000 

10,000 

178 . 29 

?l 

It 

.000 

15,000 

176.53 

If 

«! 

.000 

25,000 

170.72 

II 

It 

6,544 

ft 

.000 

50,000 

139.36 

tf 

.000 

10,000 

177.43 

ft 

.000 

15,000 

174.79 

?! 

ft 

.000 

25,000 

166.03 

ft 

« 

8,725 

.000 

50, 000 

105.63 

t! 

.000 

10,000 

176.58 

(1 

t! 

.000 

15,000 

173.05 

t* 

If 

.000 

25,000 

161.28 

I! 

ft 

21,789 

1,871.0 

50, 000 

88.70 

tf 

.000 

10, 000 

171.48 

t! 

If 

.000 

15,000 

162.64 

!! 

ft 

.000 

25,000 

129.56 

t! 

It 

14,960.0 

50,000 

91.70 

t! 

43,412 

ti 

.000 

10,000 

163.17 

1! 

.000 

15,000 

145.35 

II 

u 

12,011.0 

25, 000 

82.98 

tl 

it 

36,711.0 

50,000 

96.72 

TABLE  1  (Con't.) 


Initial 

Range 

R(0) (ft) 

Miss  Distance 
without 
control 

P0 

Miss  Distance 
with  optimal 
control 

P  (ft) 

Initial 
Veloc 1  tv 
( ft/sec*} 
v(o) 

Opt ima 1 
Steering 
Angle 
a0 

500, 000 

3,726 

.000 

10,000 

179.04 

II 

♦  t 

.  000 

15,000 

178.32 

ff 

If 

.000 

25, 000 

175.86 

If 

H 

.000 

50 , 000 

162.31 

rf 

13,088 

.000 

10, oco 

178.56 

it 

ft 

.000 

15,000 

177.43 

ft 

If 

.000 

25,000 

173-80 

rr 

tf 

.000 

50, 000 

153.01 

it 

17,450 

.000 

10,000 

178.08 

n 

ff 

.000 

15, 000 

176.65 

ft  Tl 

M 

.000 

25,000 

171.73 

If 

ff 

.000 

50, 000 

143.03 

43,573 

.000 

10, 000 

175.21 

tf 

ft 

.000 

15,000 

171.64 

ft 

.000 

25,000 

159.32 

M 

tf 

12,732.0 

50, OOC 

36 . 93 

ft 

86,  824 

.000 

10,000 

170.49 

If 

" 

.000 

15,000 

l63.il 

ff 

ff 

.000 

25,000 

138.76 

If 

ft 

56,349.0 

50,000 

91.92 
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INTERCEPT  LOGIC  AND  EQUATIONS 


To  simplify  the  bookkeeping  In  the  program  we  make 
some  changes  in  notation.  Let  x  =  x-,  ,  y  =  x_  ,  v„.  =  x0  , 

-L  ci  X  D 

j 

The  equations  of  motion  of  the  intercept  vehicle 
with  origin  at  the  target  are: 

X-,  =  X 


x2  x4 


=  a (t )  cos  a( t ) 

=  a (t )  sin  a(t ) 

where  a  (t )  -  =  ~ T 


The  thrust  angle  program,  a(t),  is  required  to 

minimize  the  time  to  intercept  (x^  =  x0  =  0)  or,  if 

that  is  impossible,  to  minimize  the  distance  to  the 

target  (min  V  2  ’ ,  2 ) .  The  application  of  Pontryagin's 

X1  +  *2 

Maximum  Principle  to  the  problem  yields  the  result  that 
a(t)  is  equal  to  some  constant  a0  .  The  equations  above 
can  then  be  integrated  to  yield: 

x^(t)  =  A(t)  cos  a°+  x^(0)t  +  x^(G) 

x2(t)  -  A ( t )  sin  a°+  x^(0)t  +  x2(0) 

x^(t)  =  A(t)  cos  a°  +  x^(0) 
x^(t)  =  A(t )  sin  a°  +  x^(Q) 
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and  denotes  the  instant  at  which  the  rocket  engine 

stops  burning. 

The  problem  is  now  one  of  determining  t°  and  o° 
such  that  x-j(t°)  =  x0(t°)  -  0  or,  if  that  is  impossible, 
such  that  -/xj(t°)  +  x|(t°)  is  minimized,  t°  being  the 
minimum  time  to  intercept  or,  if  that  is  impossible,  the 
time  to  minimum  distance  from  the  target.  If  the  coordi¬ 
nate  system  is  changed  so  that  the  intercept  vehicle  is 
initially  at  rest  at  the  origin  and  the  target  is  travel¬ 
ling  on  a  straight  line  at  constant  velocity,  then  A(t) 
represents  the  distance  of  the  intercept  vehicle  from  the 
new  origin  at  time  t.  The  equations  of  motion  of  the 
target  vehicle  in  the  new  coordinates  become: 

h  *  y3 
$2  -  y4 
y3  -  0 
=  0 


fc- 
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or 


y1:t)  =  x1(o)  +  x^(0)t 
y2(t)  -  x2(o)  +  x4(c)t 
y^(t)  =  x3(o) 

Y4  -)  =  *4(0) 

The  distance  of  th-  target  vehicle  from  the  origin  ' n 
the  new  c  'rdinate  system  is  defined  as  E(t),  or  (Fig.  5a) 

E(t)  =  y2(  t)  +  y|(t) 

which  yields  the  derivatives  to  be  used  later 
E(t)  =  j^yiU^U)  +  y2(t)y4(t)] 

1  r  "  12 

E ( t; )  =  ^[y1(t)x4(°)  -  y2(t)x3(0) 

=  x1(0)x4(0)  -  x2(0)x3(0)  ■ 

E( t)  =-  ^[y1(t)y3(t)+y2(t)y4(t)J^x1(0)x4(0)-x2(0)x3(0) J 

At  time  t  the  intercept  vehicle  can  reach  any 
point  on  the  circumference  of  a  circle  of  radius  A(t) 
by  suitable  choice  of  a.  Hence,  the  original  problem 
is  equivalent  to  one  of  determining  the  smallest  t° 
such  that  A(t°)  =  E(t°)  or,  if  that  is  impossible,  the 
t°  such  that  j A( t° )  -  E(t°)|  is  a  minimum.  The  angle  a° 
can  be  determined  from 

y-^t0)  =  E(t°)  cos  (a0  -  tt) 
y2(t°)  =  E(t°)  sin  (a0  -  v) 

Clearly,  whenever  |A(t°)  -  E(t°)J  is  a  minimum, 

[A(t°)  -  E(tc,]2  is  a  minimum  and  vice  versa.  Since 
the  latter  function  is  more  tractable  analytically,  it 
will  henceforth  be  used  to  determine  t°.  Formally  set¬ 
ting  the  derivative  of  [ A( t )  -  E(t)]2  equal  to  zero, 

^  E  A  ( t )  -  E(t)]2  3lt0  -=  2  E  A  ( t )  -  E(t)][A(tJ  -  E(  t )  ] 
yields  a  necessary  condition  for  [A(t)  -  E(t)]2  to  be  a 
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minimum;  that  Is,  A ( t )  =  E(t),  or  A(t)  E(t).  As  a 
step  in  determining  the  possible  roots  or  these  equations, 
the  quantities  t^,  tp  tp,  t^  are  defined  as  the  roots 
(if  they  exist)  of  the  equations 


or 


E(tc)  -  0 

tQ  =  "  [x.j(0)y^(0;  +  x^Ou^u'l/txqui  +  X^(O')] 

A(t1)  -  E(t1)  =  0 

A(tg!  -  B(tg)  .01  t  f 

A ( t ^ )  -  E(tp)  =  0  J  *  J 


The  quantity  t  (generally  not  equal  to  zero)  is  the 
time  at  which  the  target  vehicle  is  closest  to  the  origin 
(the  starting  point  for  the  intercept  vehicle).  The  quan¬ 
tities  tp  tp  t^  are  used  in  Chart  I  to  indicate  roots 
of  the  appropriate  equations  above,  when  they  are  known 
to  exist.  Figure  5(a)  illustrates  the  various  critical 
instants 

The  above  facts  are  used  in  Chart  I  to  determine  the 

logic  necessary  to  locate  the  solution  t°  and  to  determine 

whether  a  CATCH  (A(t°)  E(t°)]  or  a  MISS  [min[A(t°)  -E(t°)  3 2] 

is  possible.  Since  the  A  equations  change  for  t.  >  tb, 

tb  is  used  to  indicate  that  the  equations  for  t  < 

should  be  used  and  vice  versa  for  t+ .  The  times  t,,  t«, 
o  D  l  2 

t-3  and  t  may  ail  be  found  b,v  using  Newton's  method  for 

-J 

locating  the  zeros  of  a  function.  In  the  chart,  f  is 
used  to  indicate  the  function  whose  zero  is  sought.  Since 
some  of  the  functions  have  multiple  zeros  (and/or  their 
derivatives  have  zeros),  the  symbol  t  is  used  to  indicate 
a  starting  point  for  Newton's  method.  Due  to  monotone 
properties  of  the  functions  and  derivatives  involved,  use 
of  t  as  a  starting  point  guarantees  the  convergence  of 
Newton's  method. 

There  are  times  when  it  might  be  desirable  to  let 
the  intercept  vehicle  fire  at  some  instant  later  than 
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t  =  t*  (intercept  vehicle 
fires  at  this  time) 


(b)  In3  •  lal  Intercept  Coast  Trajectory 


Figure  5.  Geometrical  Configuration  of  Intercept  Trajectories 


zero,  especially  if  the  target  is  approaching  (t  >  0), 
see  Figure-  5  ( t ) ,  There  is  generally  a  maximum  permis¬ 
sible  period  of  initial  coa  t,  defined  as  t*  shown  in 
Figure  5(b).  The  remaining  material  of  this  appendix 
deals  with  the  analysis  of  this  case.  In  order  for  a 


CATCH  to  be  possible 
with  respect  to  t°. 


dt° 


A ( t °  -  t*)  =  E(t°)  . 


one  obtains 


A  o  - 
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-  t*) 


/  A-  O  \ 

n  k.  t  j 


Different iat Ing 


In  order  for  t# 
sary  conditions 


df* 

to  be  maximum  — =—  - 

dt° 

'or  the  determination 


A  ( 


o 


_  +- 


E(t( 


-  t*)  =  E(t°) 


V;  hence  two  ncces  - 
of  t°  and  t*  are 


The  problem  of  determining  t°  and  t*  can  be  simplified  i.f 

E(t)  and  E(t)  are  expressed  in  terms  of  E  ^  E(t  ), 

o  o 

minimum  distance  of  the  radius  of  target  from  th>-  int.er- 

, — - 75 - 

cept  origin,  V  =.  /x.,(0)  -t  x  j:  ( 0 ) ,  the  speed  of  the 
target  vehicle,  and  6(t),  the  angle  between  the  posi¬ 
tion  of  the  target  vehicle  at  t  and  its  position  at 
any  later  time  t,  with  6( t)  restricted  to  the  first  quadrant. 
Then 


E  /cos  0( t ) 
o'  x  ' 


E(t) 

tan  0(t)  =  V[t  -  t^] /E 
[sec2e(t)]  — ^ 


dt 


oJ 

V./E. 


o 


o 


E(t)  = 

=  V 


Eq  sin  9 ( t )  sec^@(t) 


sin  9 ( t ) 


If  At  =  t°  -  t*,  the  problem  has  become  one  of  determining 
At  and  0(t°)  such  that 


A(At)  =  E^cos  @(t°) 

A(At)  =  V  sin  0(t°) 


when  6(t°) 


is  eliminated,  the  following  equation  is 


obtained 


A  (At ) 
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it  is  easy  to  see  that,  If  A(t^)  <  V,  then  F(t)  has  at 
least  one  positive  real  zero.  If  F(w)  <  U,  then  the 
zero  must  be  less  than  t^  (P(0)  =  [VEQ]?  >0).  If 


F(t^)  >  0,  then  with  F(t)  rewritten  as 

[E- 


F(t)  =  [ v  A ( t ) ] 


A(tb)i^ 


V 


-14- 


A  ( t )  J 


t  >  t, 

—  b 


the  limit  of  F(t)  as  t  -->  «  is  clearly  negative, 
and  hence  F(t)  has  a  zero  for  t  >  t^.  If  A(tb)  >  V, 
the  Intercept  vehicle  can  fire  at  any  time  and  always 
catch  the  target  because  its  burnout  velocity  is  greaterr 
than  the  velocity  of  the  target.  Chart  II  indicates  the 
logic  steps  necessary  to  locate  the  desired  zero  of 
F(t)  . 

Once  At  is  determined,  t°  Is  specified  by 
E(t°)  ss  A(At)  with  t°  chosen  as  the  largest  root 
of  the  equation.  If  t*  =  t°  -  At  <  0,  then  intercep¬ 
tion  is  Impossible  and  the  Intercept  vehicle  should 
fire  at  t  =  0. 
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A(t2)<E(t2)  tx<t  <t, 
A(t2)>E(t2)  to<t°<t 


1 


1.2  F(t6)<0  tj>tb 

2.  F(tb)<0 

2.1  F(tb)>0  0<t1<t 

2.2  F(tg )<0  0<t1<t 


x  X  10 


Figure  6.  High  Velocity  Intercept  Paths,  Rocket  #1,  <p  =  l8l  ,  pQ  =  1745 
R(0)  =  100,000  ft. 
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Figure  7.  High  Velocity  Intercept  Paths,  Rocket  #2,  8  =  l8l  ,  pQ  =  1745  ft. 
R(0)  =  100,000  ft. 


Terminal  Error  Sensitivity  to  Initial  Measurement  and  Steering  Errors 
Rocket  #1,  V(0)  =  10,000  ft  ./sec. 


Error  Sensitivity  to  Initial  Measur 
,,  V(0)  -  25, 000  ft  ./sec. 


Terminal  Error  Sensitivity  to  Initial  Measurement 
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Figure  14.  Family  of  Interception  Paths,  Rocket  #2,  0  =  18O.50,  p  =  4363 
R(0)  =  500,000  ft.,  V(0)  =  25.000  ft. /see.. 
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Igure  15.  Terminal  Error  Sensitivity 
Rocket  #2,  V(0  =  25,000  f 


■I7I.240* 


16.  Family  of  Interception  Paths,  Rocket  #2,  <t>  =  180.5°,  po  =  4363  ft 

R(0)  =  500,000  ft.,  V(0)  =  50,000  ft ./sec. 


(FT.)  R  X 10  (FT.) 

Figure  17.  Terminal  Error  Sensitivity  to  Initial  Measurement  and  Steering  Errors, 
Rocket  #2,  V(0)  =  50,000  ft ./sec. 
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Summary  Error  Coefficient  Along  Optimal 


v<o)«2a,ooc  rt./sec. 


Figure  21.  Summary  Error  Coefficient  Along  Optimal  Path,  Hock 


Figure  23.  Steering  Angles  vs.  Tin; 
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